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SUMMARY 


This paper describes and illustrates capabilities of the FLEXSTAB Computer 
Program System. FLEXSTAB is a system of computer programs for performing 
aeroelastic analysis of a wide variety of current and future aircraft configu- 
rations. There are two versions of FLEXSTAB: a NASA control s-fixed version? 

identified as Level 1 FLEXSTAB; and an Air Force version,* ** identified as 
Level 2 FLEXSTAB, which provides for active controls analysis at low frequencies. 
The aerodynamic theory used in FLEXSTAB is applicable to both steady and 
unsteady, subsonic and supersonic flow for multiple wing-body-tail -nacelle 
configurations with a plane of symmetry. For unsteady flow calculations, an 
unsteady aerodynamic theory is used which is appropriate for the low reduced 
frequencies associated with aircraft flight dynamics. The aircraft may be 
modeled as either a rigid or flexible structure. FLEXSTAB will trim the 
aircraft in steady reference flight and compute both static and dynamic 
stability and control derivatives and the stability behavior about the trim 
condition. The airplane lifting pressure distribution, aerodynamic and inertia 
loads and deflected shape are also computed. 


INTRODUCTION 


FLEXSTAB is a system of computer programs for performing aeroelastic 
analysis of a wide variety of current and future aircraft configurations. 

There are two versions of FLEXSTAB: a NASA, control s-fixed version, identified 

as Level 1 FLEXSTAB (ref. 1 and 2); and an Air Force version, identified as 
Level 2 FLEXSTAB (ref. 3), which provides for active controls analysis at low 
frequencies. Effort is currently underway to consolidate these two versions 
into a Level 3 FLEXSTAB program. This paper describes the NASA version of 
FLEXSTAB in some detail with a brief description of the Air Force version. 


* Sponsored by NASA Ames under Contracts NAS2-5006 and NAS2-7729. 

** Sponsored by the Air Force Flight Dynamics Laboratory under Contracts 
F3361 5-72-C-l 1 72 and F-3361 5-75-C-31 32. 
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The degree to which aeroelasticity influences aircraft flight behavior, 
such as stability or controllability and structural loads, varies greatly from 
aircraft to aircraft, depending upon configuration and performance requirements. 
Obviously the aeroelastic effects become more significant with the trend to 
lighter weight, more flexible aircraft structures. The low aspect ratio, thin 
wing, slender body configurations of large supersonic aircraft are a prime 
example. The aeroelastic effects for this class of aircraft arise from complex 
structural deformation shapes (relative to higher aspect ratio configurations) 
affecting both the loads and the stability. Further, because some of their 
structural motions may have characteristic frequencies approaching those of 
their rigid body motions, dynamic stability may be significantly affected by 
the dynamics of the structure through aerodynamic coupling between the rigid 
body and structural motions. 

In recognition of these needs, NASA, and later the U.S. Air Force, sponsored 
the development of the FLEXSTAB computer programs. These are regarded as 
first generation programs which partially satisfy the need for an automated 
aeroelastic analysis tool. 

This paper presents a brief introduction to the capabilities and limit- 
ations of the FLEXSTAB programs. Some of the capabilities are illustrated by 
results from application of FLEXSTAB to actual and proposed aircraft configur- 
ations including the YF-12A aircraft, the NASA Arrow-Wing Supersonic Cruise 
Aircraft, and the Boeing Fixed-Wing SST. 


FLEXSTAB CAPABILITIES 


The NASA Level 1 FLEXSTAB is a system of fourteen separate computer 
programs (totalling about 100 000 source statements) using linear theories to 
evaluate control s-fixed static and dynamic stability, trim state, inertial and 
aerodynamic loading, and elastic deformations of aircraft configurations at 
subsonic and supersonic speeds. The functional configuration of Level 1 
FLEXSTAB is shown in figure 1. A wide range of analysis capability has been 
incorporated into the FLEXSTAB system as shown in table 1. A typical analysis 
sequence is shown in figure 2. 

The aerodynamic theory is applicable to both steady and unsteady, subsonic 
and supersonic flow for multiple wing-body-tail -nacelle configurations with 
a plane of symmetry. The aerodynamic theory is a linearized potential flow 
theory and contains a low frequency approximation to unsteady flow. A single 
aerodynamic paneling scheme can be used for steady and unsteady, subsonic and 
supersonic flow. Structural flexibility and inertial properties are modeled 
using an elastic axis representation with beam finite elements (internal to 
FLEXSTAB) .*A1 so, structural properties generated externally by programs such 
as NASTRAN^ or ATLAS can be interfaced with FLEXSTAB where a better structural 
definition is required, such as for low aspect ratio aircraft. Structural and 
aerodynamic properties are integrated using a mean-axis representation. The 
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dynamics of a flexible aircraft are resolved into structural dynamic free 
vibration modes superimposed on rigid body dynamics. The equations of motion 
are expressed for a steady, reference motion to determine trim and static 
stability, inertial and aerodynamic loading, and elastic deformations. They 
are also expressed in terms of unsteady perturbations about the reference 
motion to determine dynamic stability by characteristic roots or by time 
histories following an initial perturbation or following penetration of a 
discrete gust flow field. A brief description of these capabilities with 
corresponding limitations follows. 


FLEXSTAB AERODYNAMIC MODELING 


Geometry 

The user inputs the aerodynamic geometry of a symmetric configuration as 
a series of slender bodies with interference surfaces and thin surfaces as 
shown in figure 3. A fuselage, a nacelle, or an external store is represented 
by a slender body of revolution. The body camber is defined as camber slopes 
with respect to a straight aerodynamic mean centerline, and the boundary 
conditions are set up on a cylindrical surface of revolution. This aerodynamic 
mean centerline also defines the elastic axis for the body. Interference 
effects on a slender body are accounted for by an interference shell which is 
a cylinder of constant polygonal cross section. A wing, horizontal tail, 
vertical tail, or any other lifting surface (including control surfaces) is 
represented by linearization of the boundary conditions about a mean plane. 


Aerodynamics 

The aerodynamic influence coefficient method used in FLEXSTAB to solve 
the linearized potential flow equations is an extension of the constant 
pressure paneling technique due to Woodward (ref. 4). For each of the 
components in the aerodynamic model, FLEXSTAB computes aerodynamic influence 
coefficient matrices which relate pressure at one aerodynamic panel to the 
flow incidence and incidence rates of another panel. Incidence angles are 
expressed in terms of the configuration geometry (involving camber, twist, 
thickness, dihedral, control surface deflection), aircraft attitude and motion, 
and elastic deformations. The flow boundary conditions are the flow 
incidences. Interference effects between thin bodies and slender bodies are 
approximately accounted for by vortex panels on the thin bodies and by inter- 
ference surfaces (shells) enclosing each slender body. 

For unsteady flow calculations, an unsteady aerodynamic theory is used 
which is appropriate for the low reduced frequencies associated with aircraft 
flight dynamics. The theory is not intended to deal with high frequency 
motion such as encountered in flutter or dynamic loads due to atmospheric 
turbulence (fig. 4). The unsteady flow model is unique in that it has the 
same three-dimensional capability as the steady flow model at both supersonic 
and subsonic speeds ... a feature not found in other unsteady formulations. 
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Also, unlike other supersonic schemes such as the Mach Box method, the FLEXSTAB 
aerodynamic model does not require the introduction of "diaphragm regions" for 
wings with subsonic edges. 

Nonlinear and viscous effects are not accounted for analytically since 
the aerodynamic formulation is based on linearized potential flow theory. The 
severity of these limitations is a function of the configuration and flight 
condition. To partially alleviate this limitation, provision has been made 
for incorporating user specified changes to the aerodynamic matrices, flow 
incidence vectors, aerodynamic force and moment coefficients, and lifting 
pressure distributions. 


FLEXSTAB STRUCTURAL MODELING 


The user can model the aircraft either as a rigid or flexible structure. 
For a flexible structure, two modeling options are available (see fig. 5). If 
the structure is amenable to an elastic axis idealization, an internal 
structures program using beam type finite-elements is provided. In this case, 
the user supplies FLEXSTAB with the elastic axis geometry, and stiffness 
(El, GJ) and mass information. 

For elastic modeling of complex structures, FLEXSTAB will accept the 
structural output of external finite element programs such as NASTRAN® or 
ATLAS (fig. 6). In this case, a general three-dimensional model must be 
"reduced" to an equivalent representation consistent with the FLEXSTAB aero- 
dynamic model. That is, the flexibility and inertia definitions are reduced 
to structural nodes located on the mean aerodynamic plane of thin bodies and 
on the mean aerodynamic centerlines of slender bodies (fig. 5). 

A flexible structure can be represented in either of two ways. 

1. Static-Elastic - The structural deformations are assumed to occur 
statically, and all effects due to structural vibration are ignored. 
Using this assumption, the structural deformations are related to 
the applied loads by the flexibility matrix of the unrestrained 
vehicle and are in phase with the rigid body motion. 

2. Residual -Elastic - The structural dynamic motion is accounted for 
by using the free-vibration-mode shape amplitudes as structural 
degrees of freedom. FLEXSTAB uses the "residual flexibility" 
technique, whereby the dynamic effects of just the lower frequency 
modes (which are most likely to couple with the rigid body motion) 
are retained. The dynamic effects of the higher frequency modes are 
neglected, but their static flexibility effects are retained (hence 
the name residual flexibility). The frequencies of the retained 
modes must not violate the low frequency approximation of the 
unsteady aerodynamics. 
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In both cases the inertial relief is included in the flexibility matrix 
which is then transformed to the mean-axis system. The static-elastic approach 
is used for the trim and static stability problems; also, it is usually quite 
adequate for dynamic stability cases in which the lowest structural vibration 
frequency is relatively far removed from the frequency of the aircraft short- 
period motion. The residual flexibility technique is provided for analyzing 
dynamic stability problems in which the lower structural frequencies and the 
short-period aircraft motion may be close enough together to significantly 
couple the aircraft "rigid" body and structural dynamic motion. 

In general, the present FLEXSTAB versions are not applicable to dynamic 
loads analysis. This is because of the low reduced frequency approximation 
used for the unsteady aerodynamics and the simplifications made in the 
structural dynamic model. The structural simplifications define the fuselage 
mass and stiffness along a straight line elastic axis and neglect all rotary 
inertia terms, (fig. 7). For a wing surface, a dumbbell mass representation 
is used to represent the wing panel inertial properties. 


PROBLEM ANALYSIS 


The aerodynamic and structural representations are combined into an 
aeroelastic set of equations which govern the aircraft loads and flight 
behavior. Except for the time histories response and the postprocessing of 
the loads, solutions to these equations are computed in the Stability 
Derivatives and Static Stability (SD&SS) program. Figure 8 presents a 
schematic of the SD&SS program illustrating the user input including the 
empirical data and the program output. The basic calculations for the trim 
solution and stability and control are further described below. 


Trim Solution 

From user specified values of steady reference flight conditions (load 
factor or pitch rate, flightpath angle or thrust, yaw and roll rates, bank 
angle, and altitude or speed and dynamic pressure), FLEXSTAB computes the trim 
parameters and trimmed force coefficients shown in table 2. Three solution 
options are available: 

1. Trim solution with constant coefficients: Here, all the aerodynamic 

force coefficients are computed analytically by FLEXSTAB, and the 
trim parameters are obtained directly from the linearized equations 
of motion. 

2. Iterative trim solution: In this case the user can supply a table 

of nonlinear force and moment coefficients to replace all or a 
portion of the FLEXSTAB computed rigid aerodynamic coefficients; 
aeroelastic increments to these coefficients are computed on the 
basis of linear theory. 
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3. User specified trim parameters: This option allows the user to 

specify values for the trim parameters, thus, in general leaving 
the aircraft in an unbalanced load condition. This option is useful 
for calculating loads per unit control surface deflection, for 
specifying untrimmed flight maneuvers, or for matching rigid model 
wind tunnel test conditions. 

After the trim variables are obtained (or specified), FLEXSTAB computes 
lifting pressure distributions, aerodynamic and inertial loads at aerodynamic 
centroids, and elastic deflections due to these loads. In addition, for 
steady symmetric flight, the aerodynamic pressures can be integrated over user 
specified portions of the aircraft to obtain shear, bending, and torque 
reactions to the applied airloads (ALOADS program) when the structural model 
originates from an external finite element program. 


Stability and Control Calculations 

The FLEXSTAB analysis proceeds by perturbing the aircraft motion variables 
(both rigid and elastic) about their values for the reference flight condition. 
This results in explicit matrix equations from which the static and dynamic 
stability and control derivatives are computed. These derivatives, listed in 
tables 3 and 4, are computed for both longitudinal and lateral directional 
motions. The linearized flow equations used in FLEXSTAB govern first order 
aerodynamic effects; second-order nonlinear terms, such as the product of 
angle of attack with angle of sideslip, are neglected. For small dihedral 
configurations these second-order terms can significantly affect the yaw rate 
and sideslip stability derivatives (ref. 5). 

The stability characteristics about the reference condition are computed 
for two cases, 1. static (steady state flight), and 2. dynamic, (time varying 
motion of the aircraft). For the static case, the static stability parameters 
listed in table 2 are computed. For the dynamic case, the aircraft stability 
behavior is determined with controls fixed. (In the Air Force version of 
FLEXSTAB, reference 3, the low frequency dynamics of vehicles with feedback 
controls and sensors can be evaluated). Using the linearized equations of 
motion, the roots of the corresponding characteristic equation are computed. 
These roots in turn supply the dynamic modes of motion together with the 
associated frequencies and damping coefficients. Each real root and each 
oscillatory pair of roots are described individually. For each mode the 
following is printed: 

. Times and number of cycles to one-half (or double) and 
one-tenth amplitude 
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Frequency and period 

Logarithmic decrement and ratio of successive maximum 
displacement 


Undamped natural frequency 
Damping ratio 

Phase and amplitudes of modal coupling terms 


Time Histories Response 

The nonlinear perturbation equations of motion are numerically integrated 
using a Runge-Kutta procedure to determine the time histories response. The 
time histories response is required when the user wants to investigate the 
dynamic stability due to; (a) the large perturbations to motion variables, 
and (b) a discrete gust. Three types of discrete gust inputs may be specified, 
these are sine, 1-cosine, and modified square wave. Nonlinear aerodynamic data 
may be input as a tabular series for calculating the response. Specification 
of the discrete gust must be consistent with the limitation of the low 
frequency approximation. 


APPLICATIONS OF LEVEL 1 FLEXSTAB 


Some of the capabilities of FLEXSTAB are illustrated by the results from 
application to different airplane configurations. The NASA Arrow-Wing 
Supersonic Cruise Aircraft has been used to illustrate the prediction of 
stability and control parameters for the rigid and flexible airplane. Also, 
comparisons between FLEXSTAB results and results from other analysis tools, 
wind tunnel tests, and flight data are presented for the YF-12A, an Arrow-Wing — 
Body configuration, the Boeing Fixed-Wing SST, and the Space Shuttle Orbiter. 

The YF-12A results are for a flexible aircraft and the other three for rigid 
ai rcraft. 


NASA Arrow-Wing Supersonic Cruise Aircraft 

This configuration was used by Boeing in support of the NASA SCAR program 
(ref. 6). In that program FLEXSTAB was used to determine the aeroelastic 
loads using elastic properties supplied by the ATLAS finite-element structural 
analysis program. FLEXSTAB has also been used to predict the stability and 
control parameters for the configuration. Figure 9 shows the aerodynamic 
paneling used in the analysis. Tables 5-8 list the stability and control 
parameters for both the rigid and elastic airplane as output by FLEXSTAB. 

The elastic airplane parameters are shown for both the "free-free structure" 
and the "fixed-free structure". The free-free airplane analysis is programmed 
in the FLEXSTAB code. The elastic displacements and rotations in the free-free 
analysis are with respect to a mean axis system with the origin selected at the 
aircraft mass center. The elastic distortions relative to a mean-axis system 
correspond to those of the actual unrestrained flight vehicle and do not 
contribute to the linear and angular momenta of the relative motion with 
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respect to the body axes. Hence, the origin of the body axes remains at the 
aircraft mass center at every instant and reduces the inertial coupling 
between the overall and relative deformation motions. The fixed-free airplane 
analysis consists of using a cantilevered (clamped) flexibility matrix which 
contains implicitly the inertial relief effects. The fixed-free airplane 
analysis was patched into the FLEXSTAB code for this example to illustrate the 
significance of the mean-axis formulation. The fixed-free type of solution is 
sometimes used in present aircraft design, but the free-free solution is more 
representative of the airplane in flight. Tables 5-8 show the effects of 
using the free-free analysis. The difference between the free-free and 
fixed-free analyses for the trim a (A a =0.800°) represents the angle through 
which the mean axis rotates relative to the clamped axis. The effect of this 
rotation can be significant for the stability derivatives as presented in 
tables 6 and 7. 


YF-12A 

A comparison between the FLEXSTAB predicted and flight measured deformed 
shape of the YF-12A is shown in figure 10. These results are for a load 
factor of 1 at a Mach number of 2.8. For this flight condition the agreement 
is generally good. At subsonic Mach numbers there is not such good agreement. 
This is probably due to the chine induced vortex flows which FLEXSTAB cannot 
model . 


Arrow-Wing-Body Configuration 

Reference 7 presents comparisons of theoretical and experimental transonic 
pressure distributions for an arrow-wing— body configuration. Example results 
at M=0.85 for a = 2.1° and 7.9° are given in figure 11 for 2 spanwise 
locations. In addition to wind tunnel results, 2 sets of analytic potential 
flow results are shown; one of these is from FLEXSTAB and the other is from 
Boeing program TEA-230. The FLEXSTAB results for the individual upper and 
lower surface pressures were computed from a Boeing modified version of 
FLEXSTAB. The released version of FLEXSTAB computes only the pressure 
difference between the lower and upper surfaces (lifting pressure). 

At a = 2.1° both FLEXSTAB and TEA-230 predict pressure distributions 
which compare well with the experimental results. In contrast to the mean 
surface paneling method used in FLEXSTAB, the TEA-230 model employs 
on-the-surface paneling and boundary conditions. This feature enables the 
actual wing leading-edge pressures to be more closely predicted by TEA-230 
than by FLEXSTAB. 

At a = 7.9° the presence of a vortex spiraling off the wing leading edge 
alters the pressure distribution over a major portion of the wing. As seen 
in figure 11 the analytic results still agree with experiment at the inboard 
station but not at the outboard station. This is because neither of the two 
analytic methods can model the leading-edge vortex flow. 
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Boeing Fixed-Wing SST 


FLEXSTAB and wind tunnel comparisons for the Boeing Fixed-Wing Supersonic 
Transport (B2707-300) are shown in figure 12. These analytic results were 
computed from an early version of FLEXSTAB. The figure shows the variation 
in lift-curve slope and aerodynamic center with Mach number. Agreement 
between the measured and FLEXSTAB predicted results is good. For supersonic 
Mach numbers the aerodynamic center predicted by FLEXSTAB tends to be slightly 
aft of the measured locations, and the lift curve slope tends to be slightly 
high. 


Space Shuttle Orbiter 

Figures 13 and 14 show comparisons of pitching moment coefficients due to 
pitch rate c m g (table 3) and due to rate of change of angle of attack 
c m £ (table 4) q for the Space Shuttle Orbiter. Wind tunnel values for the 
sum (cmq + c m £) are compared to FLEXSTAB computed values in figure 13. The 
agreement is remarkably good considering the unstreamlined shape of the 
Orbiter. In contrast to wind tunnel measurements in which only the sum 
c mq' + c m $ is generally measured, FLEXSTAB computes each derivative separately. 
These are shown in figure 14 where it is seen that c m £ is destabilizing at 
the supersonic Mach number. 


LEVEL 1 FLEXSTAB DOCUMENTATION, AVAILABILITY AND IMPLEMENTATION 


The documentation for Level 1 FLEXSTAB is listed under references 1 and 2. 

Both a CDC and an IBM version of the NASA FLEXSTAB Program are available 
from COSMIC through a lease arrangement. Inquiries concerning the lease 
should be made to: 

Computer Software Management and Information Center 
Barrow Hall 
University of Georgia 
Athens, Georgia 30601 

The CDC version has been executed on the following computers, using 
FORTRAN Extended compilers: 

1. CYBER 70, Model 73, SCOPE 3.4 Operating System 

2. CDC 6600, KRONOS 2.1 and SCOPE 3.4 Operating Systems 

3. CDC 7600, SCOPE 2.1.3 Operating System 
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The IBM version requires both the G and H compilers and has been executed 
on the following computers: 

1. IBM 360-91/OS 

2. IBM 360-67/OS 

Additional information on the program implementation requirements is 
summarized in reference 1. 


LEVEL 2 FLEXSTAB CAPABILITIES 


The Air Force has developed additional capabilities for FLEXSTAB for the 
evaluation of the low frequency dynamics of control configured vehicles (ref. 
3). These capabilities are subject to the same limitations as Level 1 FLEXSTAB 
imposed by the structural modeling and the low frequency aerodynamic approxi- 
mation. The additional capabilities can be summarized as: 

. Active Controls Capability 

User Specification of Linear Control System 
. Sensor Simulations 

. Time Histories of Controlled Airplane 
. RMS Aircraft Motion Due to Turbulence PSD 
. Frequency Response Calculations 

. Structural Damping 

. Modal Truncation Formulation for the Structural Dynamic Model 
. Plotting of Dynamic Analysis Data 
. Data Interface For Flight Simulators 
. Fore and Aft Gust Capability 

. Aerodynamic Hinge Moments 

This extended capability of FLEXSTAB allows the analytical representation 
of arbitrary combinations of active control systems and linear systems analysis 
of the response to turbulence and control commands. Nonlinear response to 
discrete gusts and control commands may also be evaluated by time histories 
generated by the system. 

The functional configuration of the Level 2 FLEXSTAB system is shown in 
figure 15. The program system totals twelve separate programs with about 
140 000 source statements. Five of the programs remain the same as for Level 1 
FLEXSTAB, six have been modified, mainly to allow definition of motions at 
the sensor locations used by the flight control system, and one program, the 
linear system analysis, is new. 

The program has been used by the Air Force for research on the following 
aircraft; YF-16, B-l, C-5A, and the F-lll TACT. Also the program is being used 
in a preliminary design mode of operation for system integration in support 
of fire control systems on the airborne laser weapons program. 
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Figure 16 shows a typical result from Level 2 FLEXSTAB for the B-52E air- 
craft response to a lateral gust with and without a yaw stability augmentation 
system (yaw SAS). 


LEVEL 2 FLEXSTAB DOCUMENTATION, AVAILABILITY AND IMPLEMENTATION 


The documentation for Level 2 FLEXSTAB is listed under reference 3. The 
program is available only for application to Department of Defense problems 
from the Air Force by sending a request to: 

AFFDL 

Flight Control Division 
Wright-Patterson Air Force Base 
Ohio 45433 

The Level 2 version has been executed on the CYBER 74 and CDC 6600 under 
the KR0N0S 2.0 and 2.1, SCOPE 3.4 and NOS/BE Operating Systems. 


CONCLUDING REMARKS 


The development of the Level 1 FLEXSTAB computer program system provides 
a first generation program to partially satisfy the requirement for an auto- 
mated aeroelastic analysis tool. The system is based on integrated aero- 
dynamic, structural, and dynamic analytical methods applicable to aircraft 
configurations having a plane of symmetry. The aerodynamic analysis may be 
supplemented or modified with empirical data. Also the system includes the 
low frequency aerodynamic effects appropriate for evaluating the stability of 
large aircraft. 

Through application to several different aircraft configurations it has 
been shown that the linearized theory employed is capable of providing aero- 
elastic analysis of complex configurations. The linearized theory provides 
good correlation with experiment as long as the flow over the configuration 
remains attached, although in some cases it may also provide useful data for 
separated flow conditions. Further work is required to develop methods of 
representing leading-edge vortex flow. 

The Level 2 system integrates the aerodynamic and structural methods used 
in Level 1 with control system dynamics. The resulting system is capable of 
analyzing the static and low frequency dynamic characteristics of control 
configured vehicles. 

By refining the structural modeling and providing full unsteady aero- 
dynamic capability, the FLEXSTAB system could be applied to prediction of 
dynamic loads at higher structural frequencies. 
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Symbols used in Table 1 and Figure 8 


El 

beam flexural rigidity 

GJ 

beam torsional rigidity 

[LSC] 

steady lifting aerodynamic influence coefficient matrix 

P, Q, R 

components of angular velocity, co 

U, V, w 

perturbation components of velocity 

p 

density of air 

6 , 

Euler angles 


Subscripts: 

1 reference state 

Superscripts: 

• (dot) first derivative with respect to time 
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TABLE 1.- PROGRAMS OF THE LEVEL 1 FLEXSTAB SYSTEM 






Data Flow 

No. 

Program Title 

Program 

Name 

Purpose 

Input 

Output 




AIRPLANE MODELING SECTION 


1 

Geometry 

definition 

program 

GD 

To geometrically construct an airplane 
configuration that complies with the 
requirements of the aerodynamic and 
structural theories 

• Basic geometric description 
of the airplane (wing plan- 
form, airfoil shapes, wing 
twist, body radius distribu- 
tion, body camber, etc.) in 
its cruise or unload condition 

• Geometric defining quantities 
required for a linearized aero- 
dynamic and structural repre- 
sentation (such as panel corner 
points and location of line 
doublet segments) 

CalComp plots of the airplane's 
configuration 

2 

Aerodynamic 

influence 

coefficient 

program 

AIC 

To produce matrices that relate the 
change in pressure on a panel to the 
change in flow incidence at another 
panel 

A- = 

and to produce matrices that relate 
the change in pressure on a panel to 
the rate of change in flow incidence 
at another panel 

a p i 

• Mach number 

• Geometric defining quantities 
produced by the GD program 

• Aerodynamic influence coefficient 
(AIC) matrices 

3 

AIC matric 

correction 

program 

CAIC 

To correct the [LSC] AIC matrix 
using 

(a) a correction matrix [C2] based on 

wind tunnel or flight-test data 

(b) a replacement technique in which 
the user specifies elements of the 

matrix 

• The [LSC] AIC matrix 

• A corrected [LSC] matrix 


TABLE 1.- CONTINUED 


\Jo. 

Program Title 

Program 

Name 

Purpose 

Data Flow 

Input 

Output 

STRUCTURAL MODELING SECTION 


4 

Internal structural 
influence 

coefficient 

program 

ISIC 

To produce flexibility matrices based 
on beam theory and consistent with 
paneling used for computing the AIC 
matrices 

< 

i 

•Geometric defining quantities 
produced by the GD program 
•Elastic axis geometric 
description 

•Els and GJs of beam segments 
•Mass distribution 

•Flexibility matrices (Static- 
elastic) 

•Mass matrix 

5 

Elastic axis plot 
program 

EAPLOT 

To plot the structural elastic axes 

•Elastic axis data 

•CalComp plots of the 
airplane's elastic axis 

6 

Normal modes 
program 

NM 

To calculate the free-vibration 
normal mode shapes and then to 
construct normal mode matrices 

•Number of normal modes 
requested 

•Structural (ISIC) flexibility 
matrices 

•Mode shape matrices 
•Residual flexibility matrices 

7 

Normal modes plot 
program 

NMPLOT 

To plot the normal mode shapes 

•Normal mode shapes 

•CalComp plots of the mode 
shapes 

8 

External structural 
influence 
coefficient 
program 

ESIC 

To accept a flexibility matrix from 
a source external to the FLEXSTAB 
system and by interpolation produce 
a new flexibility matrix consistent 
with paneling used for computing 
the AIC matrices 

•Geometric defining quantities 
produced by the GD program 
•Structural mode shapes 
•External flexibility matrix 

•New flexibility matrices 
•Mass matrix 
•Normal mode matrices 
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TABLE 1.- CONCLUDED 


No. 

Program Title 

Program 

Name 

Purpose 

Data Flow 

Input 

Output 

PROBLEM ANALYSIS SECTION 

9 

Stability derivatives 
and static stability 
program 

SD&SS 

To calculate the static and 
dynamic stability 
characteristics of the 
airplane 

•Geometric defining quantities 
produced by the GD program 

•AIC matrices 

• Flexibility matrices 
•Normal mode matrices 

• Flight condition (altitude, p, 

R-p Q*|, 

•Wind tunnel data 

• Exterior downwash data 

•Trim parameters (a*, 6 e 
T 1 ,/3 1 ,6 ar and5 ri + 

•Deformed shape 

•Steady and unsteady pressure 
distribution 

•Static and dynamic stability 
derivatives 

•Static stability parameters 
(static margin, neutral point, etc.) 

•Control effectiveness 

10 

Pressure distribution 
plot program 

PDPLOT 

To plot the steady and 
unsteady pressure 
distributions 

•Steady or unsteady pressure 
distribution 

•CalComp plots of the pressure 
distribution 

11 

Time histories 
program 

TH 

To help the user evaluate 
the stability of an airplane 
using large-disturbance 
theory (nonlinear equations) 

•Geometric data 

• Dynamic stability derivatives 

• Stability data (o^, /J-j, 

7i, ete .) 

•Wind tunnel data 

•Values of motion variables 
versus time of flight (u, v, w, 
u, V, w, 0,0, i//, 0,0. ^ etc.) 

12 

Time histories plot 
program 

THPLOT 

To plot the time histories 
(motion variables versus 
time) of the flight 

• Motion variables versus time 

•CalComp plots of the time 
histories 

13 

Structural loads 
program 

S LOADS 

To print the structural 
elastic axis loads 

• Structural load matrices 

•Structural loads on elastic 
axis 

14 

Aerodynamic loads 

ALOADS 

To calculate air loads 

• Air load matrices 

•Air loads 


TABLE 2.- STATIC STABILITY AND TRIM PARAMETERS 


TRIM PARAMETERS 

a 

angle of attack 

a 

angle of sideslip 

5 e 

elevator deflection angle 

5 a 

aileron deflection angle 


rudder deflection angle 

T 

thrust 

7 

flightpath angle 

TRIMMED FORCE COEFFICIENTS 

c L 

lift coefficient 

C D 

drag coefficient 

C m 

pitching moment coefficient 

C y 

sideforce coefficient 

c c 

rolling moment coefficient 

C n 

yawing moment coefficient 

STATIC STABILITY PARAMETERS 

h n h 

static margin 

h n 

neutral point 

h m h 

maneuver margin 

a 6e «n(=-i a ) 

elevator angle per g (turn, pullup) 

35 e /3V 

speed stick stability 


K> 

O' 

a\ 


TABLE 3.- STATIC STABILITY DERIVATIVES AND CONTROL EFFECTIVENESS COEFFICIENTS 



lift coefficient at a = 5 e = 0 
drag coefficient at a = 5 e = 0 
pitching moment coefficient at a = 6 e = 0 
lift coefficient due to angle of attack 

drag coefficient due to angle of attack 

pitching moment coefficient due to angle of attack 

lift coefficient due to pitch rate 
drag coefficient due to pitch rate 
pitching moment coefficient due to pitch rate 
lift coefficient due to elevator deflection 

drag coefficient due to elevator deflection 

pitching moment coefficient due to elevator deflection 

sideforce coefficient due to roll rate 

rolling moment coefficient due to roll rate 



yawing moment coefficient due to roll rate 
sideforce coefficient due to yaw rate 
rolling moment coefficient due to yaw rate 
yawing moment coefficient due to yaw rate 
sideforce coefficient due to sideslip 

rolling moment coefficient due to sideslip 

yawing moment coefficient due to sideslip 

sideforce coefficient due to aileron deflection 
(antisymmetrically deflected) 

rolling moment coefficient due to aileron deflection 
(antisymmetrically deflected) 

yawing moment coefficient due to aileron deflection 
(antisymmetrically deflected) 

sideforce coefficient due to rudder deflection 
rolling moment coefficient due to rudder deflection 
yawing moment coefficient due to rudder deflection 




TABLE 4.- DYNAMIC STABILITY DERIVATIVES AND CONTROL EFFECTIVENESS COEFFICIENTS 


'L A 

L a 


c d£ 

OL 

C m A 


a 


O. A 

L U 

C DJ 

C K 

Cl A 

q 

C °q 

C^A 

m q 

C,,A 


CpA 

P 

%* 

c vf 

C n^ 


lift coefficient due to angle-of-attack rate 
drag coefficient due to angle-of-attack rate 

pitching moment coefficient due to angle-of-attack rate 
lift coefficient due to speed change 
drag coefficient due to speed change 
pitching moment coefficient due to speed change 
lift coefficient due to pitch acceleration 
drag coefficient due to pitch acceleration 
pitching moment coefficient due to pitch acceleration 
sideforce coefficient due to roll acceleration 
rolling moment coefficient due to roll acceleration 
yawing moment coefficient due to roll acceleration 
sideforce coefficient due to yaw acceleration 
rolling moment coefficient due to yaw acceleration 
yawing moment coefficient due to yaw acceleration 


c 4 

CpA 

V 

C n A 

n P 



sideforce coefficient due to sideslip rate 

rolling moment coefficient due to sideslip rate 

yawing moment coefficient due to sideslip rate 

lift coefficient due to elevator deflection 
rate 


Cq • drag coefficient due to elevator deflection 
"e rate 

C m - pitching moment coefficient due to elevator deflection 
°e rate 


•Cy • sideforce coefficient due to aileron deflection 
"a rate (antisymmetric deflection) 

Cp. rolling moment coefficient due to aileron deflection 
"a rate (antisymmetric deflection) 

C n . yawing moment coefficient due to aileron deflection 
"a rate (antisymmetric deflection) 



sideforce coefficient due to rudder deflection 
rate 


Cp ■ rolling moment coefficient due to rudder deflection 
°r rate 

C n . yawing moment coefficient due to rudder deflection 
6 r rate 


( 

i 


N5 

O' 


TABLE 5.- NASA ARROW-WING SUPERSONIC CRUISE AIRCRAFT STATIC STABILITY 

AND TRIM PARAMETERS 



Free-Free Airplane 

Fixed-Free Airplane 


Rigid 

Static-Elastic 

Static-Elastic 

Trim parameters 




deg. 

5.318 

5.441 

6.241 

(3-| deg. 

0 

0 

0 

6 ei deg. 

2.298 

14.156 

14.156 

6 ai deg. 

0 

0 

0 

6 fl deg. 

0 

0 

0 

lb. 

144 000 

144 000 

144 000 

7*| deg. 

7.319 

7.349 

6.549 

Static stability parameters 




h n -h 

0.0983 

0.0183 

0.0183 

h n 

0.6213 

0.5413 

0.5413 

h m' h 

0.1010 

0.0205 

0.0205 

36 e /3n(=^ p ) deg./g 

-8.000 

•2.416 

-2.421 

35 e /3v deg./ft/sec. 

0.0070 

0.0010 

0.0010 


Weight = 298 636 Kg 
(657 000 1b) 


C L = 0.0666 C.G = 0.523 c 

M = 2.9 h = 17 435 m 

(57 200 ft) 
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TABLE 6.- NASA ARROW-WING SUPERSONIC CRUISE AIRCRAFT STATIC STABILITY 
DERIVATIVES AND CONTROL EFFECTIVENESS COEFFICIENTS 



Free-Free Airplane 

Fixed-Free Airplane 


Rigid 

Static-Elastic 

Static-Elastic 

C Lq 

-0.090394 

-0.087729 

-0.093639 

C D 0 

0.004326 

0.004246 

0.004306 

C m 0 

0.016662 

0.010725 

0.010831 

C. /deg. 

0.029008 

0.026645 

0.023918 

Cn /deg. 

0.003038 

0.002916 

0.002955 

C /deg. 
" Q 

-0.002851 

-0.000486 

-0.000438 

Ci A /rad. 
L n 

0.636205 ' 

0.345434 

0.500099 

H 

CrjA/rad. 

a 

0.068107 

0.047567 

0.072209 

^m A ^ r ^d • 

m q 

-0.881881 

-0.727361 

-0.730188 

Ci /deg. 

0.000939 

0.000568 

0.000804 

(j 0 

CD^/deg. 

0.000059 

0.000046 

0.000780 

C ^, e /deg - 

-0.000918 

-0.000566 

-0.000570 

C^A/rad. 
y n 

-0.002286 

-0.001445 

0.006297 

H 

Ci A/rad. 

v n 

-0.133523 

-0.099685 

-0.100626 

H 

C n A/rad. 

P 

0.006054 

0.005165 

0.004863 

C V A/rad. 

y r 

0.219438 

0.128576 

0.172671 

C^/rad. 

-0.007389 

-0.018225 

-0.025462 

C-Vrad. 

"r 

-0.167888 

-0.114434 

-0.122142 

Cy^/deg. 

C^/deg. 

C n i3 /de 9- 

-0.003305 

-0.001722 

-0.002574 

0.000038 

0.000214 

0.000327 

0.001279 

0.000292 

0.000443 

c v 5a /deg - 

-0.0 

0.000010 

0.000013 

C„ /deg. 

-0.000115 

-0.000057 

-0.000057 

C„ 5 ' a / d «9- 

0.000011 

0.000003 

0.000004 

C V 6r /deg - 

0.000228 

0.000071 

0.000163 

S/ de 9- 

-0.000002 

0.000015 

0.000003 

C n /deg. 
o r 

-0.000172 

-0.000069 

-0.000086 
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TABLE 7.- NASA ARROW-WING SUPERSONIC CRUISE AIRCRAFT DYNAMIC 

STABILITY DERIVATIVES 



Free- Free Airplane 

: ixed-Free Airplane 


Rigid 

Static-Elastic 

Static-Elastic 

C|_*/rad 

-0.048826 

0.064902 

0.244706 

CrWrad 

-0.006662 

0.004098 

0.024072 

C m^ /rad 

0.058374 

-0.021734 

-0.024995 

C L A 

0.064077 

0.068511 

0.045412 

C Dy 

0.005153 

0.059110 

0.004428 

C m A 

,n u 

0.008657 

0.007856 

0.008273 

C LA /rad 

-0.280092 

0.345434 

0.500099 

o 

Q> 

CL 

-0.015575 

0.047567 

0.072209 

C m£ /rad 

-0.062581 

-0.727361 

-0.730188 

o 

< 

Q) 

CL 

-0.021108 

-0.011722 

-0.020624 

Cg^/rad 

0.009710 

-0.003486 

-0.002221 

C n ^/ ra d 

0.014685 

0.008969 

0.010627 

C yf/ rad 

-0.043381 

-0.031054 

-0.040741 

Cj^/rad 

-0.043000 

-0.009016 

-0.008496 

c nf/ rad 

-0.043619 

-0.051572 

-0.049689 

CyA/rad 

-0.059523 

-0.058859 

-0.059655 

P 

Cg^/rad 

-0.002412 

-0.001590 

-0.001705 

C n^ rad 

-0.015210 

-0.015543 

-0.015375 


TABLE 8.- NASA ARROW-WING SUPERSONIC CRUISE AIRCRAFT LONGITUDINAL AND LATERAL DIRECTIONAL 

DYNAMIC CHARACTERISTICS 


Mf = 2.9, h! =17 435 m, 7i = 7.3°, eg = 0.523c, wt = 298 636 Kg 
(57 200 ft) (657 000 lb) 

V = 1015 km/hr 



/ 

Short Period 

Phugoid 



Period 

(sec) 

T 1/2 or T 2* 
(sec) 

C i /2 Of C 2 * 
(cycles) 

Period 

(sec) 

T 1/2 or T 2* 
(sec) 

C 1 /2 or C 2 * 
(cycles) 

Free-Free 

Rigid 

3.55 

Ti/ 2 - 2.44 

C -j 12 = 0.687 

368.3 

T 2 = 28 904 

C 2 = 78.5 

Airplane 

Static- 

elastic 

8.54 

T-|/2 = 2.69 

C-j /2 = 0.316 

214.9 

T-,/2 = 267.8 

C-j/2 = 1-25 

Fixed-Free 

Airplane 

Static- 

elastic 

8.99 

Ti/ 2 = 2.84 

C-jy2 = 0.316 

215.9 

T 1/2 = 225.7 

C-j/2 = 




Dutch roll 

Roil 

convergence 

Spiral 

mode 



Period 

(sec) 

T 1/2 orT 2* 

(sec) 

C i /2 or C 2 * 
(cycles) 

\m 

Phase 

angle 

(deg.) 

T 1/2 
(sec) 

t 1/2 or T 2 
(sec) 

Free- Free 

Rigid 

4.83 

t 1/2= 10 - 7 

C 1/2 = 2.21 

0.22 

41.45 

T 1/2 = 1.36 

T 2 = 308.0 

Airplane 

Static- 

elastic 

12.15 

t 1/2 = 241. 3 

C-j/2 = 18.85 

5.56 

-144.43 

T 1/2 = 1.49 

T 2 = 60.7 

Fixed- Free 
Airplane 

Static- 

elastic 

10.23 

T 2 = 107.93 

c 2 =10.5 

6.17 

-146.74 

t 1/2 = 139 

T 2 = 61.2 


*Time or cycles to half or twice amplitude. 


Job 

control 


Problem 

specifications 


Host 

Operating 

System 


Interfaced Systems 

• NASTRAN* 

• ATLAS 

• SAMECS 


[ Structural j 

( Loads j 

l data * — | 

1 \ data *- 


AIC matrix correction 


Aerodynamic influence 
coefficients 


Geometry definition 


Aerodynamic Modeling 


External structural 
influence coefficients 


Normal modes plots 


Normal modes 


Elastic axis plots 


Internal structural 
influence coefficients 


Structural Modeling 


Aerodynamic loads 


Structural loads 


Time histories plots 


Time histories 


Pressure distribution 
plots 


Stability derivatives 
& static stability 


Problem Analysis 


Figure 1.- Level 1 FLEXSTAB functional configuration. 



Figure 2.- FLEXSTAB analysis sequence 




X Axis 


Thin surfaces 

• Wing, tails, struts 
Bodies of revolution 

• Fuselage, nacelles 
Interference shell 

(Polygon of constant 
cross section) 



COMPONENTS ASSEMBLY 



Figure 3.- FLEXSTAB aerodynamic geometry modeling. 


Increasing complexity of math model 


Rigid Static Dynamic 

Flexibility Flexibility 



Figure 4.- Limitations on low frequency aerodynamic approximation. 
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• INTERNALLY PRODUCED ELASTIC AXIS MODEL 



Straight line elastic axis 
for fuselage 


• EXTERNALLY SUPPLIED FINITE ELEMENT MODEL REDUCED TO LINE AND PLANAR MODEL 



Translational degrees 



Figure 5.- FLEXSTAB structural modeling. 


FLEXSTAB 


ATLAS, NASTRAN • 


• Aerodynamic model and 
equations of motion for 
unrestrained elastic 
aircraft 


• ^cj - Free flexibility 

• Aeroelastic loads 




centroids; reactions 


• Structural model 





c l - Clamped flexibility 
'5| - Mass distribution 


H 

N - 1 


I * J “ Generalized stiffness 
mass and mode shape 
matrices 



Figure 6.- FLEXSTAB/ ATLAS, NASTRAN® 


interface. 
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TRUE AERODYNAMIC GEOMETRY 



Upswept aft fuselage 


FLEXSTAB AERODYNAMIC GEOMETRY 



Fin placed at edge 
of interference shell 


Fuselage is a body of revolution plus 
camber slopes defined along straight 
aerodynamic mean centerline. 


FLEXSTAB ELASTIC AXIS 


Mass and stiffness properties 
Defined along aerodynamic 
mean centerline. 

(No rotary inertia) 


Figure 7.- FLEXSTAB modeling limitations. 



/ SD&SS data deck 

• Altitude 

• Rotational velocities; P-j, Q-j, R-j 

• Load factor n,, instead of Q* 

Z 1 1 

• Roll angle; 0^ 

• Wing geometry; S w , c, b 

• Flightpath angle, 7^, or thrust, T -j 

• Rigid, static-elastic, or 
residual-elastic 

• Type of reference motion 

• Number of engines and location 

• Engine angular momentum 

• Control surface data 

• Panel area ratios 

• External pressure distribution 

• Recycle data 



Airplane trim parameters 



Figure 8 


Schematic of the SD&SS program 





Figure 9.- NASA Arrow-Wing supersonic cruise aircraft aerodynamic peneling. 


: 3 === 


Nacelle (£ 
”Y WBL 100 
Body £ 




Figure 10.- FLEXSTAB/ flight deformed shape of YF-12A. 
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Figure 11.- Flat wing pressure distribution of Arrow-WLng — body. 
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Figure 12.- Effect of Mach number on lift-curve slope and 
aerodynamic center of Boeing fixed-wing SST. 



Figure 13.- Space shuttle orbiter values of C ~ + C^. 
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Figure 14.- FLEXSTAB values of and Cm^* 



Figure 15.- Level 2 FLEXSTAB functional configuration 



Side acceleration Yaw Rate rad/sec Rudder deg 

m/sec^ (+Nose Rt.) 


159 090 Kg 648.6 Km/Hr 1219.2 m 

(350 000 lbs) (350 kn EAS, M = 0.569) (4000 ft) 



Figure 16.- B-52 response to lateral gust at B.S.616. 
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